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I.   Introduction 

It  has  long  been  knoi^m  that  certain  flows  of  liquids 
with  free  boundary  siirfaces  are  unstable.   For  steady  flows 
this  means  that  a  small  perturbation  introduced  at  a  point 
of  the  flovj  produces  a  dist\.irbance  which  increases  in  the 
flovj  direction  until  it  completely  alters  the  notion.   This 
?nechanism  iijas  used  by  Rayleigh  to  describe  the  breakup  of  a 
steady  .iet  of  liquid  into  drops.   In  the  case  of  unsteady 
floi-js,  instability  implies  that  certain  initial  rserturbations 
will  increase  in  time  until  they  significantly  change  the 
flow.   Taylor  showed  that  the  plane  free  boundary  of  a  semi- 
infinite  liquid  is  unstable  if  the  liquid  moves  v;ith  constant 
acceleration  normal  to  the  surface,  with  the  acceleration 
directed  into  the  liquid.   Previcusl;-  Penney  and  Pryce  had 
shovm  that  the  surface  of  a  pulsatin.c:  underwater  explosion 
bubble  is  like^^7ise  unstable.   In  this  case,  also,  the  fluid 
moves  with  accelerated  motion  normal  to  the  surface. 

With  these  results  in  mind,  we  wish  to  account  for  the 
formation  of  drops  resulting  from  the  breakup  of  accelerated 
plane,  cylindrical  and  spherical  layers  of  liquid.   These 
phenomena  occur  whenever  exulosives  arc  detonated  In  or  near 
liquids.   T?or  example,  the  spray  dome  produced  by  an  under- 
water explosion  results  from  a  layer  of  vrator  throv;n  into 
the  air  by  the  shock  \\ra.ve    from  the  explosion.   The  present 
theory  attemiots  to  ds scribe  the  breakup  of  this  layer  into 
drops  of  spray.   As  a  secoid  example,  consider  the  production 
of  an  aerosol  (i.e.  dispersion   of  liquid  drops  In  air)  by  a 
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bomb  containinp:  an  explosive  surrounded  by  a  layer  of  liquid. 
The  high  pressure  gases  resulting  frora  the  explosion  propel 
the  liquid  layer  outward,  and  it  Liltinately  breaks  ut)  into 
drops.   The  present  theory  also  attempts  to  account  for  this 
brealcup  . 

In  our  proposed  m.echanlsm  of  breakup  v;e  assuiTie  that 
there  is  a  zero  order  accelerated  flow  of  the  liqiiid  laj^er. 
The  first  order  perturbation  of  this  flow,  which  satisfies 
linear  equations,  is  represented  as  a  series  or  integral  of 
normal  modes.   Some  of  the  modes  are  found  to  bo  unstable 
in  the  sense  that  their  a'apli tildes  increase  indefinitely 
with  time.   Among  these,  one  or  more  are  m.ost  unstable,  in 
the  sense  that  they  grow  most  rapidly.   In  the  case  of 
constant  acceleration  the  unstable  modes  grow  exponentially, 
and  the  most  unstable  modes  have  the  largest  exponents. 
The  most  unstable  modes  are  assumed  to  be  responsible  for 
the  breakup  of  the  laycv    into  pieces,  which  then  become 
spherical  due  to  surface  tension.   The  number  of  pieces 
depends  upon  the  number  of  positive  and  negative  regions 
into  vjhich  the  most  unstable  mode  function  is  decom.poscd 
by  its  nodal  lines.   Thus  it  is  possible  to  estimate  the 
number  and  size  of  the;  resulting  drops. 

In  the  follovjing  section  (II)  the  problem  is  formulated 
for  the  case  of  a  plane  layer.   In  section  III  the  zero 
order  solution  is  considered,  and  in  section  IV  the  first 
order  perturbation  is  obtained.   This  perturbation  solution 
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is  analyzed  in  section  V,  where  the  formation  cf  drops  Is 
discussed.   Practical  conclusions  ctnd  rccora'tiendations  are 
presented  in  section  VI.   Corresponding  results  for 
cylindrical  and  spherical  layers  ^^Jill  be  given  in  another 
paper. 

II.   Foi-'r!'iulation 

¥e  consider  the  inotion  cf  a  plane  slab  or  layer  of 
liquid  which  separates  tv;o  gases  at  different  pressures. 
The  liquid  villi  be  assumed  to  be  incoripressible  and 
inviscid,  and  the  motion  of  the  gases  vjill  be  neglected. 
The  initial  positions  'ind  shapes  of  the  bounding  surfaces 
viill  be  prescribed,  as  xvell  as  the  initial  velocity  of  the 
liquid.   The  problem  is  to  determine  the  subsequent  positions 
and  shapes  of  these  surfaces,  and  the  motion  of  the  liquid. 

Assuming"  that  the  motion  is  irijational,  the  velocity 
u(x,y,z,t)  of  the  liquid  i«  derivable  from  a.  potential 
<f;(X,y,2,t)  vjhich  satisfies  Laplace's  equation.   Thus 


u  =  -  V<^ 


2.  s7^^, 


The  pressure  p(x,y,z,t)  in  the  liquid  can  be  expressed  in 
terms  of  ^,    the  constant  density. p,  and  the  acceleration  of 
.rrravity  g  by  Bernoulli's  equation 


pp"   ==  (tj^  -  gz  -  ^  [V^f 


It  is  assumed  that  the  positive  z  axis  points  vortically 
upward . 
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The  boundary  surfaces  are  assumoo  to  be  given  by 

h-  z  =  F._(x,y,t)        i  =  1,2  . 

The  subscript  i  =  1  corresponds  to  the  urper  s\.ir.face.   The 
gas  pressiire  on  both  surfaces  is  assuoied  to  be  independent 
of  X  and  y,  and  is  denoted  by  p.  (i  =  1,2).   Then  the 
dynamic  and  kinematic  boundar;''  conditions  ma^"  be  written 
respectively  as 


on  z  =  F. 


(i  =  1,2) 
^-    *z-  *x^x-  Vly^  ^"it  =  ^    on  z  =  F.       (i  =  1,2)   . 

I'i"'-  (^),  T.  is  the  surface  tension  of  the  surface  i.   The 
initial  conditions  will  be  introduced  when  they  are  needed. 
If  tne  flow  is  bounded  by  a  cylindrical  tube  parallel 

to  the  z-axis,  the  condition  r—  =  0  is  imposed  on  this  tube. 

on 

If  the  lower  surfa,ce  z  =  F^  is  a  rigid  surface  with 
proscribed  motion,  the  above  formulation  applies  provided 
(<)  is  omitted  for  i  =  2.   TherX^re  this  problem  will  also 
be  considered,  since  it  is  so  simi\ir  to  the  original  problem. 

s 

III.   Zero  Order  Solution 

If  the  bounding  surfaces  are  initially  planes  normal 
to  the  z-axis,  and  if  the  initial  velocity  is  parallel  to 
the  z-axis  and  independent  of  x  ixnCi   y,  the  solution  will  be 
independent  of  x  :..r\d.   y.   Then,  if  u(t)  denotes  the  z-comnoncnt 
of  velocity,  wc  have  from  (1,  2) 


( 
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(f)'^   =   -   u(t)z   +   b(t) 


nO     o. 


Since    F.''   depends    on  t    only,   vjc   liavo   from   (5-7)  »    denoting  t 
differentiation   by   a   dot, 

8.  b   =    (u+  g)P° +i  u^+  p'^p^.  i   =   1,2 

9.  ^  =  ^it  '-   "   -'2 

Frora  (9)  we  have,  if  ;i  denotes  the  initial  separation 
between  the  surfaces 

10.  F2  =  F^  -  h    . 

From  (8),  eliminating  b,  we  obtain  using  (10) 

11 .  (u+  g)ph  =  Pp  -  p^ 

Equation  (11)  j'-ields  u,  and  therefore  u  if  its  initial  valiie 
is  given,  provided  Pp  -  p-,  is  constant  or  a  known  function 
of  t  (or  even  of  u  and  t).   If  the  pressure  difference  is 
constant,  u  Is  constant,   Hox-jevcr  if  Pp  -  P-.  depends  upo^n  Fp, 
then  (11)  becomes  a  second  order  equation  for  Fp(t). 
Finally  b  is  determined  from  (8),  except  for  an  inessential 
additive  constant. 

In  the  problem  in  which  Fp(t)  is  a  prescribed  moving 
rigid  surface,  u  and  F-,  are  given  in  terms  of  Fp  by  (9,  10), 
and  b  by  (5)  i;itfi  i  =  1,  while  (5)  for  i  =  2  is  omitted, 
as  well  as  its  consequence  (11). 
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IV.   First  Order  Perturbation 

Suppose  that  the  initial  boundary  surfaces  differ 
slightly  from  planes,  and   or  that  the  initial  velocity 
differs  slightly  from  the  constant  z-Component  assumed  in 
the  above  solution.   Then  we  may  expect  the  subsequent 
solution  to  differ  slightly  from  the  above  solution,  which 
we  will  henceforth  call  the  zero  order  or  unperturbed 
solution.   If  e  is  a  measure  of  the  maximum  deviations  of 
the  initial  data,  we  assume  that  the  solution  may  be  written 
in  the  form 


12 .       <|)  =  <t)°  +  e  (j)-^  +  .  .  .  ,     F.  =  F?  +  e  F^ 


+ 


From  (2,  12)  we  see  that  ^      also  satisfies  Laplace's 
equation.   Nov;  inserting  (12)  into  (5,  6)  we  obtain  for  the 
first  order  perturbation  4'  and  F   the  equations 

13.  ^^-{n+g)F^+u^^-  {-l)^p'^1^^^F^=  0     onz  =  F?    (i=l,2) 

II+.  ct)^+F^  =  0   on  z=  F?    (i=l,2)  . 

In  (13,  ll-i )  and  the  sequel  the  superscript  one  on  <j)  and  F 
is  omitted. 

To  solve  Laplace's  equation  for  <})  we  assume  that  (f)  is 
a  product  of  a  function  of  z  and  t  multiplied  by  a  function 
'HXjy)  .   We  find 

15.  ^  =  [g^{t)e^^  +  g2(t)e-^^]^l;(x,y) 

16.  {^^  +    k^)\\/   =   0 
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If  the  flovj  is  unbounded  in  the  x,y  plane,  then  (l6)  has  a 
bounded  solution  for  every  real  value  of  k,  and  the  general 
solution  for  <}>  is  an  integral  with  respect  to  k  of  the 
expression  on  the  rioiht  of  (15),  in  ^^^hich  g,  and  gp  depend 
upon  k. 

However,  if  the  flow  is  bounded  by  a  rigid  cylinder 
parallel  to  the  z-axis,  and  if  this  cylinder  intersects  the 
x,y  plane  in  a  curve  C,  then  the  function  vji  must  also  satisfy 
the  condition 

17.  |l  "  ^     °^^  ^    • 

If  C  is  closed,  then  (16,  1?)  have  solutions  only  for  a 

2  2 

denumerable  set  of  positive  values  of  k  ,  and  for  k   =  0. 

In  this  case  the  general  solution  for  <})  is  a  series  of  terras 
of  the  form  given  by  (1?),  summed  over  the  various  values 
of  k. 

In  order  to  complete  the  solution  wo  must  find  g^ ,  gp , 
P,  and  Fp .   To  this  end  we  assume 

18-  P.  =  f^(t)>^(x,y).  i  =  1,2   . 

Now,    inserting   (15,    l3)    into    (1-^,    lL|.)    and  m.aking  use    of    (16) 
we    obtain 

.      kP'5^      ^      -kF°  kP?  -kP? 

19.      g-j_e         +g2®         "  -  (u+g)  f^  +  ukg^e      -^-ukgpO        •■- 


+  (-1)  Vp-^T.f.    =   C         (i    =   1,2) 
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kF?       -kF? 
10.  kg^e   ^  -  kg^^         ^  +  f^  =  0  (i  =  l,2) 


■'^rom  (20)  we  find  g-,  and  gp  in  tcrnivS  of  f^^  and  f  o : 

21.  g,^   =   -k-^e    l(f^-e--^  f2 )  (1  -  e-^^''^") "" 

22.  g2  =  k-^e   l(f^-e^h  f^)  { 1  -  e2kh)-l 


Novr  using    (21,    22)    in    (19),    and   combining   the    ti-jo   equation 
( 19) ,   we    obtain 


23.         (f^-e-k^   f2)+k(u-.g)(f,-e-kl^f2) 


+   k'V^o    -(f3_+  T-(_^T2e  f  2  ^    "    ° 


2IL.         (V   -e^h   f:,)-k(u+p)(f,  -e^h   f^) 


-   k-\^p"^(f^ +T^-^T2e^^^^   f^) 


Equations  (23,  2I4.)  are  a  pair  of  linear  horaogenecus 
second  order  ordinary  differential  equations  for  f-,  and  f p . 
The  coefficients  are  constant  only  if  u  is  constant,  which 
case  we  shall  treat  in  detail.   First  however,  we  wish  to 
Tjoint  out  the  limiting  forms  which  (23,  2I4.)  attain  whtn  h 
becomes  infinite,  namely 

25.  *f\  +  k(u+g+ k-T-^p"-^)fT_  =  C 


26.  £'2   -   k{n+g  -k^T2P~'^)i^2   ^   ^ 
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Either  of  these  is  the  equation  for  the  perturbation  of  the 
surface  of  a  semi-infinite  region  of  liauid,  in  agreement 
with  Taylor's  result.   V/e  note  that  for  sufficiently  large 
constant  u  one  of  the  above  equations  will  have  a  growing 
exponential  solution  while  the  other  will  not. 

If  u  is  constant  we  seek  solutions  of  (23,  21^.)    of  the 
form 

27.  f^  =  A^e^^  i  =  1,2   , 

Inserting  (27)  into  (23,  2J|)  and  simplifying  leads  to  an 
equation  for  a   and  another  for  A^/A-,  ,  namely 


23.    a  +  a'"k^p"-'-(T-,  +  To)coth  kh 


-  k^(u+g+ k^T^p"-'-)(u+g-  k^Tpp"-^)  =  0 


29.  A^/A^  =e^^'^(a'+  k[u+P:]  +  k^T^o"-^)(a^  + k[u+£T]  -k^Tgp"^)""'-   . 

The  solutions  of  (23)  for  a"  are  riven  by  (v;hen  T-,  =  Tp) 

30.  ph^TT-'-c^  =  -x^  coth  X  t   x(p^  +  x'-^-  sinli'^  x)"^/^ 


The  abbreviations  used  in  (30)  are 


31.        X  =  kh   ,     G  =  ph^TJ-'-lu+gl 


Since  there  are  four  values  of  a,  each  f.  is  a  si.im  of 
four  terms  of  the  form  in  (27). 

In  the  problem  in  which  P'pv^)  ^^  ^^^^  prescribed  motion 
of  a  plane  rigid  surface,  F2(t)  is  zero.   All  equations 
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throup-h  (2/4.)  apply,  but  (13,  l5,  19)  hold  only  for  i  =  1  and 
fp  =  0  in  (21,  22).   Instead  of  (23,  2i+)  we  obtain  for  f^ 
the  equ'ition 

23'.        *f^  +  k  tanh  kh  (u+g  +  k^T-j_p"")  f-j_  =  0 

For  infin5.te  h  this  also  has  the  limiting  form  (25). 

V.   Instability  and  Drop  Formation 

2 

If  a   is  negative,  a   is  pure  imaginary  and  the 

corresponding  ex^3onential  is  oscillatory.   This  occurs  v/hen 

the  negative  sign  is  chosen  in  (30).   Only  when  the  plus 

2 

sign  is  chosen  is  there  a  range  of  x  for  which  g   is  positive, 

and  thus  one  value  of  c  is  real  and  positive  leading  to  a 

growing  exponential.   In  this  range  there  is  a  greatest 

value  of  a  corresponding  to  a  particular  value  of  x  which 

we  v;ill  call  x    ,  since  the  exponential  co7'-resPonding  to 
max'  '  '      ^ 

this  value  of  x  grows  most  rapidly. 

To  find  X    we  first  note  that 
max 

-lo  Ij-   .  ,  -2      T  .. 

32.  X   sinii    X  <  1.3 


Thus  if  6  »  JT73   =   1.1^,  (30)  becomes  (using  the  plus  sign) 

33-  ph-'^T^^-'-a'^  =  x(p-x^  ccth  x) 

2  1/2 

From  (33)>  g  is  positive  for  x  between  zero  and  P  '  , 

approximately.   It  has  its  maximijim  at 
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p  1/2 

3k'  X    =  {4)  =  h 

■^^  max    ^3 


p  lu+p:  i   ^ 
3T, 


Rocallins;  that  by  (31),  x  =  kh,  we  may  introduce  k    which 

max 

is  the  value  of  k  corrcspondin'T  to  x    ,  and  thup  to  maximum 

"-■  max ' 

instability.      From   (3^-)    wo   have 


35.  k^       =  %t£i.       . 

-^'^  max  3T^ 

To  find  the  number  of  pieces  into  which  the  layer 

breaks  up  due  to  the  growth  of  the  term  corresponding  to 

k    ,  vjt'  must  examine  the  function  ^I/(x,y)  corresDondinc:  to 
max '  -TK     ,j  /  .      ^ 

k    .   It  is  clear  that  if  the  layer  is  imbounded  in  the 
max 

x,y  plane  there  will  be  many  functions  \['{k,j)    corresponding 

to  k    .   On  the  other  hand,  if  the  region  is  bounded  (1?) 

applies  and  there  may  be  no  solution  t  corresponding  exactly 

to  k    ,  although  there  will  be  m.any  functions  4/  corres- 

ponding  to  values  of  k  near  k    .   Thus  in  both  cases  there 
^  max 

will  be  many  functions  -^  with  exactly  or  approximately  the 

samxC  value  of  k  {7:^    k    )  and  therefore  all  the  corresponding 

terms  will  grow  at  about  the  same  (maximum)  rate. 

Consequently  the  exact  manner  of  breakup  will  depend  upon 

the  extent  to  v;hich  ther-c  various  terms  are  excited  by  the 

initial  perturbation. 

Nevertheless  it  is  possible  to  make  a  rough  estimate 

of  the  mamber  of  pieces  into  vjhich  the  layer  breaks.   First, 

in  the  case  of  an  unbounded  region,  one  solution  \|/ 

correspond inr-  to  k     is 

max 
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36. 

where 

37. 


\1/  =  0 


2T:i{x/\T_+  y/x^) 


Wc  now  consider  the  regions  of  the  x,3^  plane  in  which 
the  real  part  cf  \}/  is  positive  (or  negative).   These  are 
rectangles  bounded  by  nodal  lines,  and  their  dimensions 

are  -^  and  -^ ,  with  area  — r— ,   If  \,  and  \p  are  related 

t>y  (37)  J  the  area  of  such  a  rectangle  is  minimized  when 

2r/? 
^n  =  \2  ~  k *   'T^'^^'-S  the  minimum  area  of  such  a  rectangle  is 


max 


33 


^1 


2-^' 


max 


,  2 

DTi;  T- 


o  1  u+g 


It  is  to  be  expected  that  (38)  will  represent  roughly  the 
area  of  a  piece  into  which  the  lay^r  breaks  \xp    since  all 
parts  of  the  surface  in  this  region  move  in  the  same 
direction.   The  volume  of  a  piece  will  be  given  by  (3S) 
multiplied  by  the  thickness  h.   The  radius  r  of  the  sphere 
into  which  the  piece  will  ultimately  deform  is  then 
(us Ing  ( 11) ) 


39. 


r  = 


grcT^h 


1V^ 


'-2piu+g!  J 


o^T^h^  ^1/3 


L2iP2-p^ 


The  same  results  are  obtained  if  the  layer  is  bounded  oy   a 
tube  of  rectangular  cross-section,  and  they  will  also  apply 
for  tubes  of  almost  any  other  shane . 
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VI.   Conclusions  and  Recommendations 

On  the  basis  of  (39)  three  means  are  available  for 
diminishing  drop  size,  namely 

a)  Reducin,^  the  layer  thicknet's  h 

b)  Reducing  the  surface  tension  1\    by  using  wetting 
agents 

c)  Increasing  the  explosive  pressure  Pp . 

Of  those  the  first  two  seem  reasonable,  while  the  last 
appears  Impractical. 

VJe  note  from  (39)  that  diminishing  h  has  a  twofold 
effect:   the  acceleration  is  greater  and  the  pieces  are 
thinner.   It  vjcxild  seem  that  the  reduction  in  drop  size, 
due  to  decreasing  h,  with  the  consequent  reduction  in 
settling  out  of  the  dror>s  would  mcrcj  than  compensate  for 
the  diminution  of  liquid  in  the  bomb,  for  a  certain  range 
of  h.   The  theory  seems  to  bo  too  crude  to  allow  calculrition 
of  an  optimum  h,  but  this  could  be    explored  experimentally. 

The  incorporation  of  wotting  agents  in  the  bomb,  either 
mixed  with  the  liquid  or  in  a  thin  layer  separated  by  a 
fragile  membrane,  seems  promising.   If  only  one  surface  can 
be  wet  it  should  be  the  inner  surface. 

If  the  formula  (39)  is  i\sed  to  calculate  the  size  of 
drops  produced  by  a  bomb,  it  is  necessary  to  use  a  value  of 
h  less  than  the  ori^rlnal  value  duo  to  radial  expansion  of 
the  liqv.id  layer,  and  its  consequent  tninning.   One  might 
guess  that  the  layer  radius  increases  tenfold,  and  the 
thickness  h  consequently  diminishes  100  fold  before  breakup 
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occurs.   More  precise  resLilts  will  be  /riven  in  the  forth- 
coming paper  en  the  breakup  of  spherical  and  cylindrical 
layers . 

To  estimate  the  time  required  for  breakup  to  occur, 
wo  assume  that  the  surface  oorturbation,  given  essentially 
by  (27),  roaches  the  value  h.   This  gives  for  the  breakup 
time 


ho. 


t  = 


max 


-1    n 
log  J 


I  ' '^  13 


1/ii 


.OS" 


h 
A 


log  J 


m  (1+0)  A  is  the  initial  perturbation  of  the  surface,  a 


max 


was  obtained  from  (3.3,  3l|),  and  the  last  equality  follows 
from  (11) . 
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